Measurement of finite-frequency current statistics in a single-electron transistor 
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Electron transport in nano-scale structures is 
strongly influenced by the Coulomb interaction 
which gives rise to correlations in the stream 
of charges and leaves clear fingerprints in the 
fluctuations of the electrical current. A complete 
understanding of the underlying physical pro- 
cesses requires measurements of the electrical 
fluctuations on all time and frequency scales, 
but experiments have so far been restricted to 
fixed frequency ranges as broadband detection 
of current fluctuations is an inherently difficult 
experimental procedure. Here we demonstrate 
that the electrical fluctuations in a single electron 
transistor (SET) can be accurately measured on 
all relevant frequencies using a nearby quantum 
point contact for on-chip real-time detection of 
the current pulses in the SET. We have directly 
measured the frequency-dependent current 
statistics and hereby fully characterized the 
fundamental tunneling processes in the SET. 
Our experiment paves the way for future inves- 
tigations of interaction and coherence induced 
correlation effects in quantum transport. 

The electrical fluctuations in a nano-scale conduc- 
tor reveal a wealth of information about the physical 
processes inside the device compared to what is avail- 
able from a conductance measurement alone. ^ ^ Zero- 
frequency noise measurements are now routinely per- 
formed using standard techniques, but more recently the 
detection of higher-order current correlation functions, 
or cumulants, has attracted considerable attention in 
experimentaPti^ and theoretical^M^ studies of charge 
transport in man-made sub- micron structures. Measure- 
ments have for example been encouraged by the intrigu- 
ing connections between current fluctuations and entan- 
glement entropy in solid-state systems^^ as well as by 
the possibility to test fluctuation theoremP^JI^ at the 
nano-scale which are now a topic at the forefront of non- 
equilibrium statistical physics. 

Experiments on current fluctuations have mainly fo- 
cused on zero-frequency current correlations^ while 
measurements of finite-frequency cumulants of the cur- 
rent have remained an outstanding experimental chal- 
lenge. However, noise detection in a restricted frequency 
band gives only partial information about correlations 
and measurements that cover the full frequency range 
are necessary to access all relevant time scales that char- 



acterize the transport process. In a zero-frequency mea- 
surement, correlation effects are integrated over a long 
period of time and important information about charac- 
teristic time scales are lost. To observe the dynamical fea- 
tures of the correlations, finite- frequency measurements 
are required.^ ^ 28 34 j^^^g even been shown theoretically 
that the frequency-dependent third cumulant (the skew- 
ness) of the current may contain additional information 
about correlations and internal time scales of a system 
compared to the finite- frequency noise alone, for instance 
in chaotic cavitie^^ and diffusive conductors."^ 

Figure [l^-c shows a typical time trace of the currents 
in our nano-scale SET consisting of a quantum dot (QD) 
coupled via tunneling barriers to source and drain elec- 
trodes. The QD is operated close to a charge degeneracy 
point in the Coulomb blockade regime, where only a sin- 
gle additional electron at a time can enter from the source 
electrode and leave via the drain electrode. The applied 
voltage bias eV is larger than the electronic temperature 
ksT which prevents electrons from tunneling in the op- 
posite direction of the mean current. While a stream 
of electrons is driven through the QD, a separate cur- 
rent is passed through the nearby quantum point contact 
(QPC) whose conductance is highly sensitive to the pres- 
ence of additional electrons {/S.q = 0, 1) on the QD.^^^ 
By monitoring the switches of the current through the 
QPC we can thus detect in real-time single electrons tun- 
neling through the left (L) and right {R) barriers and 
thereby determine the corresponding pulse currents II (t) 
and Init). The time trace of the currents illustrates the 
electrical fluctuations of interest here. 

We now analyze the current fluctuations obtained 
from data measured over 24 hours during which 300 
million electrons passed through the SET. Our high- 
quality measurements allow us to develop a complete 
picture of the noise properties of the SET, which does 
not only focus on the zero-frequency components of the 
fluctuations, but contains the full frequency-resolved 
information about the noise and higher-order corre- 
lation functions. From the measured finite- frequency 
noise spectrum we extract the correlation time of the 
current fluctuations. The noise spectrum, however, is 
a one-frequency quantity only which does not reflect 
correlations between different spectral components of 
the current. To observe such correlations, we employ 
bispectral analysis and consider the finite-frequency 
skewness (or bispectrum) of the current. The skewness 
shows that the current fluctuations are non-gaussian 
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FIG. 1: Single electron transistor (SET) and finite- frequency noise spectra. a,b, Fluctuations of the current through 
the left (red curve) and right (blue curve) tunneling barriers of the SET. c, Atomic force microscope topography of the SET 
consisting of a Coulomb blockade quantum dot (QD) coupled via tunneling barriers to left (L) and right (K) electrodes. The 
scale bar corresponds to 250 nm. A bias difference between the electrodes drives a stream of electrons through the QD from 
L to R. An electron entering the QD from L causes a suppression of the current through the nearby quantum point contact 
(QPC) until the electron tunnels out of the QD to R. The occupation of the QD (A^ = 0, 1) is inferred from the time-dependent 
current through the QPC (shown with black in a). The corresponding pulse currents through the left (red) and right (blue) 
barriers are discretized in time steps of At = 40 /is and shown displaced for clarity, a. In b, the same currents are shown on a 
much longer time scale and discretized in time steps of At (grey spikes) and 20 At (full lines), respectively d. Noise spectrum 
(blue curve) of the current through the right barrier (i?i?) and cross-correlations (green curve) between the two currents (Li?). 
The thickness of the lines indicate experimental error estimates. Model calculations (dashed lines) of the spectra are in excellent 
agreement with the experiment using F^ = 13.23 kHz and Vr = 4.81 kHz for the tunneling rates across the barriers and a 



detector rate oiVn 



= 0.3 MHz. The mean dwell time of electrons on the QD is Td — F^ :^ 210 /xs and the correlation time is 
55 /is. The inset shows a schematic model of the charge transport process. 



on all relevant time and frequency scales due to the 
non-equilibrium conditions imposed by the applied 
voltage bias. Our measurements are supported by model 
calculations which are in excellent agreement with the 
experimental data. The results presented here provide a 
fundamental understanding of the electrical fluctuations 
in SETs which are expected to constitute the basic 
building blocks of future nano-scale electronics. 

Results 

Noise spectrum. The measured noise spectrum for 
the right tunneling barrier is presented in Fig. [T[i. The 
finite-frequency current-correlation function S)^^{ijS) is 
defined in terms of the noise power a^ 

where /a(co'), a = L^R^ is the Fourier transformed cur- 
rent and double brackets ((...)) denote cumulant aver- 
aging over many experimental realizations. Figure [l]i 
shows the Fano factor F^'^\uj) = S^l^{uj) / el for the right 
tunneling barrier with the mean current / = {Ihif)) = 



{Init)) being constant in the stationary state. The noise 
is symmetric in frequency, S^^j^{uj) = 6'^^(— cj), and re- 
sults are shown for positive frequencies only For uncor- 
related transport the noise spectrum would be white, i. 
e., F'^'^\uj) = 1 on all frequencies, corresponding to a 
Poisson process. Our measurements, in contrast, show a 
clear suppression of fluctuations below the Poisson value 
at low frequencies. This is due to the strong Coulomb in- 
teractions on the QD which introduce correlations in the 
stream of electrons propagating through the SET: each 
electron spends a finite time on the QD during which it 
blocks the next electron entering the QD. The measured 
noise spectrum has a Lorentzian shape whose width is 
determined by the inverse dynamical time scale of the 
correlations. Our measurements of the finite-frequency 
noise thereby enable a direct observation of the correla- 
tion time Tc — 55 /j.s of the transport, Fig. [T]i. 



To corroborate our experimental results we calculate 
the noise spectrum of the schematic model in the inset 
of Fig. [l]i. Single electrons tunnel from the left electrode 
onto the QD at rate Tl and leave it via the right electrode 



3 



10- 



N 
X 





n 

\ experiment / 
Ik / 

I ■ ■ ■ ■ 


y theory >i 

♦ 

* — . . 


^^^^^^^^^^^^^^^^ 



1.0 
0.8 
P3)0.6 
0.4 
0.2 



-10 



F(3) 



\ @ \ LRR 

-®"'--.,\ i©"; 



RRR 




0.4 



0.2 



0.0 



■■■■< © 



\":-..©i 

\ ■■■■© •! 



-10 -5 5 
0)7271 (kHz) 



10 



an 



F<2)(0) 




-5 5 
(0/271 (kHz) 



10 



-10 -5 5 
(0/271 (kHz) 



10 




-5 5 
(0/271 (kHz) 



FIG. 2: Finite- frequency skewness. a, Experimental results and model calculations of the third Fano factor (the frequency- 
dependent skewness) F'^^\(j0i^U2) — S^^\(j0i,uj2)/e'^I- Experimental results are shown above the diagonal cji = (jJ2 and model 
calculations below. Contour lines indicate F^^"^ = 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8. The parameters for the model calculations are 
given in the caption of Fig. [l] b, Symmetries of the bispectrum S^^\uji,uj2). Several important symmetry conditions follow 
from the definition, including symmetry with respect to interchange of frequencies S^^\(jJi,(jJ2) = S^^\(jU2,u;i). Knowledge of 
the bispectrum in any of the regions (T) — © is sufficient for a complete description of the bispectrum. We have measured the 
skewness in region (2) and the cross-bispectrum (Fig. [3| in regions (T) — (3). c,d,e. Finite- frequency skewness along the three 
blue lines in a. Full lines are experimental results, while dashed lines show model calculations. The thickness of the full lines 
indicate the experimental error estimates. The third order Fano factor F^^\(jJi,(jJ2) approaches the zero-frequency shot noise 
F^'^\0) for (jj2 = and IcJil being much larger than the inverse correlation time r""*^ as seen in c. 



at rate Tr. The Fano factor is theiPl 



where theoretically Tc = {Tl + T^^)"^ is identified as the 
correlation time. This expression qualitatively explains 
the measured noise spectrum. Quantitative agreement 
is obtained by also taking into account the finite detec- 
tion rate Tjj of the QPC charge sensing scheme. The 
three parameters Tl^ Tr^ and Tjj can be independently 



extracted from the distri bution o f waiting times between 
detected tunneling events. The model calculations 
(see Methods) are in excellent agreement with measure- 
ments over the full range of frequencies, demonstrating 
the high quality of our experimental data. 

Cross-correlations. Cross-correlations between 
the left /l(^) and the right lR{t) currents can also 
be measured, Fig. [Tli. For the schematic model, 

— (2) 

the (cross-correlation) Fano factor reads F^^{uj) = 

I + r%)/{[rL + Tr]^ + Lu^), which 



Re[S^^l{cv)]/eI = (r 
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FIG. 3: Finite-frequency cross-bispectrum. a, Measurements of the frequency-dependent cross-bispectrum ^^^^(cji, C(;2). 
The cross-bispectrum shows a reduced symmetry compared to the bispectrum and measurements in the regions ©-®, Fig.|^, 
are required for a complete characterization of the cross-bispectrum. For comparison, the bispectrum F^^\(jO 1,0)2) from FiglSa 
is shown as a semi-transparent surface above the cross-bispectrum. b, Bispectrum and cross-bispectrum along the line c<;2 = 0, 
where they coincide, c, Bispectrum (full blue line) and cross-bispectrum (dashed red line) along the line cji = 0. 



in the zero-frequency limit coincides with the noise spec- 
trum, i^^^(O) = F'^^^(O), as a consequence of charge con- 
servation on the QD. The two currents are clearly cor- 
related at frequencies that are lower than the inverse 
correlation time r~^, but the cross-correlator eventu- 
ally reaches zero at higher frequencies. Interestingly, 
the cross-correlations of the detected pulse currents are 
slightly negative at high frequencies. This is due to the fi- 
nite resolution of the QPC charge sensing protocol which 
is not able to distinguish current pulses that are sepa- 
rated in time by an interval which is shorter than the 
inverse detector rate F^^. 

Higher-order cumulants. We now turn to measure- 
ments of higher-order finite-frequency cumulants. The 
m'th finite-frequency current correlator corresponding to 
a time-dependent current I{t) is defined as 

((/(cJi) • • • i{uJm))) = 27r(5(a;i + . . .^UJm)S^'^\uJu . . . ,CJ^-l) 

where translational invariance in time implies frequency 
conservation as indicated by the Dirac delta function S{uj) 
and S^'^\uJi, . . . ^(jJm-i) is the polyspectrum.^^'^ In the 
case m = 2, the polyspectrum yields the noise power 
spectrum S^'^\uj)^ while the skewness (or bispectrum) 
is given by m = 3 with the corresponding Fano factor 
F(3)(cJi,cj2) = S^^\uJi,uj2)/e^L We focus here on the 
frequency-dependent skewness 5^^^ (cJi, (^2), although our 
experimental data in principle allows us also to extract 
cumulants of even higher orders. 

The measured skewness. Fig. [2^, shows a much richer 
structure and frequency dependence compared to the 
noise spectrum. The skewness obeys several symme- 



tries following from the definition. Fig. [SJd. We exploit 
the mirror symmetry with respect to interchange of fre- 
quencies, S^^\uJi^uj2) = S^^\uj2j'^i)j to compare mea- 
surement and model calculations: experimental results 
are presented above the diagonal ooi = while model 
calculations are shown below. Fig. [2^. The nonzero 
bispectrum indicates non-Gaussian statistics on all fre- 
quencies and shows strong correlations between different 
spectral components of the current. Importantly, these 
correlations are not a consequence of non-linearities in 
the detection scheme, but are solely due to the physical 
non-equilibrium conditions imposed by the applied volt- 
age bias. The pulse currents are directly derived from 
the tunneling events and the influence of external noise 
sources, including the amplification of the QPC current, 
is thereby explicitly avoided. Intuitively, one would ex- 
pect the correlations to vanish, if the observation fre- 
quency is larger than the average frequency of the trans- 
port. Surprisingly, however, a certain degree of corre- 
lation persists even if one frequency is large, while the 
other is kept finite. This is in stark contrast to the sec- 
ond Fano factor F^'^\uj) which approaches unity in the 
high-frequency limit, Fig. Id, corresponding to uncorre- 
lated tunneling events. 

For the schematic model in Fig.[l]i, the finite- frequency 
skewness readJ22 



J 



^n-=i[(rL + r^P + a;f]' 

having defined 7^ = Tl ^ T% 7I = 3{Tl^Tr)^, and 
^3 = co'i + C(J2. This expression qualitatively explains 



5 



the measured finite-frequency skewness and siiows tiiat 
tiie skewness iias a complex structure wiiicii is not just 
a simple Lorentzian-shaped function of frequencies. We 
note that the third Fano factor F^^^ (cji, co'2) reduces to 
the zero-frequency limit of the noise F^^^(O) in the limit 
002 = and \uji\ T~^. This is immediately visible 
in Fig. |2]3. Again, quantitative agreement is achieved 
by including the finite detection rate F^:) in the model 
calculations as illustrated explicitly in Figs. ^-e. The 
analytic expression shows that the frequency dependence 
of the skewness, unlike the auto- and cross-correlation 
noise spectra, is not only governed by the correlation 
time Tc = {Vl + T^^)"^. The skewness has an involved 
frequency dependence, given by several frequency scales, 
which cannot be deduced from the noise spectrum alone. 

Higher-order cross-correlations. In Fig. 

[3^, we finally consider the cross-bispectrum 

^Li?i?(^i'^2) = Re[5'[^^^(a;i,cj2)]/e^/, measuring 
the third-order correlations between tunneling electrons 
entering and leaving the SET. The three non-redundant 
permutations of the left and right pulse currents lead 
to a reduced symmetry of the cross-bispectrum, as 
visualized in comparison with the auto-correlation 
bispectrum in Fig. |3^. The cross-bispectrum coincides 
with the auto-bispectrum on the line CJ2 = 0, where 
they approach the zero-frequency correlation of the shot 
noise for \uji\ > r~^. Fig. pp. In contrast, for \uj2\ 7^ 0, 
the cross-bispectrum shows a frequency dependence 
similar to the cross-correlation shot noise and the 
anti-correlations due to the detection process become 
visible. Fig. 

Discussion 

We have measured the current statistics of charge 
transport in an SET and directly determined the dynam- 
ical features and time scales of the current fiuctuations. 
From the measured frequency-dependent noise power we 
extracted the correlation time of the fiuctuations. The 
noise power is a single-frequency quantity only and in 
order to investigate the correlations between different 
spectral components of the current using bispectral 
analysis we measured the frequency-dependent third 
order correlation function (the skewness). The skewness 
shows that the current statistics are non-gaussian on all 
frequencies due to the applied voltage bias. Our experi- 
mental results are supported by model calculations that 
are in excellent agreement with measurements. The re- 
sults presented here are important for future applications 
of SETs in nano-scale electrical circuits operating with 
single electrons. Our accurate and stable experiment 
also facilitates several promising directions for basic 
research on nano-scale quantum devices. These include 
experimental investigations of fiuctuation relations at 
finite frequencies and higher-order noise detection of 
interactioiP^ and coherencd^ induced correlation effects 
in quantum transport under non-equilibrium conditions. 



Methods 

Device fabrication. The device was fabricated by local 
anodic oxidation techniques using an atomic force micro- 
scope on the surface of a GaAs/AlGaAs heterostructure 
with electron density 4.6 x 10^^ m~^ and a mobility of 64 
m^/Vs. The two-dimensional electron gas residing 34 nm 
below the heterostructure surface is depleted underneath 
the oxidized lines, allowing us to define the quantum dot 
(QD) and the quantum point contact (QPC). 
Measurements. The experiment was carried out in a 
^He cryostat at 500 mK with an applied bias of 900 /iV 
across the QD in order to ensure unidirectional transport 
and to avoid the influence of thermal fluctuations. The 
QPC detector was tuned to the edge of the flrst conduc- 
tion step. The current through the QPC was measured 
with a sampling frequency of 500 kHz. The tunneling 
events were extracted from the QPC current using a step 
detection algorithm and converted into time-dependent 
pulse currents: Time was discretized in steps of At = 40 
/is and in each step the number of tunneling events An in 
(out) of the QD was recorded. The current into (out of) 
the QD at a given time step is then I{t) = {—e)An/ At. 
Error estimates. We estimated the errors of the mea- 
sured spectra by dividing the experimental data into 30 
separate batches. The spectra were determined for each 
batch individually and their standard deviation was used 
as a measure of the experimental accuracy. 
Finite- frequency cumulants. To deflne the 
flnite-frequency cumulants of the current we con- 
sider the m-time probability distributiorP^^ 
p(^)(72i, ti; . . . ; n^, t^) that Uk electrons have been 
transferred during the time span [0, t/c] for /c = 1, . . . , m. 
The corresponding cumulant generating function is 

(x, t) = log I E ^^"^ ^)^"''' I 

with X = (Xi,---,Xm), t (ti,...,t^), and n 
(ni, . . . , rim)' (An equivalent deflnition uses only a single 
but time-dependent counting fleld.^^) The m-time cumu- 
lants of P^^\n]i) are then 

((n(ti) • • • n(tra))) = di^, ■ ■ ■ aix„^(™)(x,*)lx^o, 

with the corresponding m-time cumulants of the current 

((/(ti) • • • J(U))) =dt,--- dt^ {{nih) ■ • • n(t^))). 

These equations deflne the cumulant averages denoted 
by double brackets ((...)). In the Fourier domain, the 
current cumulants can be expressed as 

((/(C^l) • • • IM)) = 2^(5(c^i+. . .^CUm)S^'^\cJu . . . ,C^m-l) 

since the sum of frequencies is zero in the stationary state. 
The Fourier transformed current is denoted as I{uj) and 
S^'^\uJijUJ2j ' ' ' ,^m-i) is the polyspectrum,^^ which for 
m = 2 and m = 3 yields the noise spectrum S'''^\uj) 
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(second cumulant) and the bispectrum (thhd cumulant) 
S''^\uJi^uj2) respectively. We note that the noise spec- 
trum S^'^\uj) in the stationary state is a single-frequency 
quantity which can be related directly to the one- 
time probability distribution P^^\n^t) via MacDonald's 
formuld^a 5(2) (cj) = uj dt sin {ujt)f^{{r?)){t). The bis- 
pectrum S^^^ (co'i, 6^2) in contrast is a two-frequency quan- 
tity which reflects correlations of the current beyond 
what is captured by P^^\n^t) alone. 
Theoretical model. We consider the probability vector 
\p{^)) = [Poo(^),Pio(^),Poi(^),Pii(^)]^, where the first in- 
dex denotes the number of (additional) electrons on the 
QD, i = 0, 1, and the second index denotes the detected 
number of (additional) electrons on the QD as inferred 
from the current through the QPC, j = 0, 1. The proba- 
bility vector evolves according to the rate equation 

^\p{xL,XR'^t)) = M(xL,Xi?)|p(xL,Xi?;^)), 

having introduced separate counting fields xl and xr 
that couple to the number of detected electrons that have 
passed the left and the ri ght bar riers, respectively.— The 
matrix M{xl^Xr) reads ^SmUl 
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where factors of e*^^(^) in the off-diagonal elements cor- 
respond to processes that increase by one the number 
of detected electrons that have tunneled across the left 
(right) barrier, see Fig. [l]3,d. The detector rate of 
the QPC detector scheme tends to infinite for an ideal 
detector only, but is finite otherwise. 
Calculations. For the calculations of the finite- 
frequency cumulants it is useful to write the matrix as 

M(XL, xr) = Mo + (e^^^ - 1)Il + (e^^^ - 1)Ir, 



where Mq = M(0, 0) and lL(i?) is the super operator for 
the detected current through the left (right) barrier.^^ 
Additionally, we need the stationary state |0)), found by 
solving Mo|0)) = and normalized such that ((0|0)) = 1, 
where ((0| = (1, 1, 1, 1). We then define the projectors^ 
P = |0))((0| and Q = 1 - P and the frequency- 
dependent pseudoinverse^'^ Il{uj) = Q[iuj -\- Mo]~"^Q, 
which is well-defined even in the zero-frequency limit 
^ 0, since the inversion is performed only in the 
subspace spanned by Q, where Mq is regular. These 
objects constitute the essential building blocks for 
our calculations. The frequency-dependent second 
and third cumulants, S'''^\uj) and ^^^^ (c^i, 6^2), can then 
be evaluated following Refs. |47|48j and [33], respectively. 
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